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We present a theory of magnetotransport through a system of two coupled electronic orbitals, 
where the electron spin interacts with a (large) local magnetic moment via an exchange interaction. 
For the physical realization of such a setup we have in mind, for example, semiconductor quantum 
dots coupled to an ensemble of nuclear spins in the host material or molecular orbitals coupled to a 
local magnetic moment. Using a semiclassical approximation, we derive a set of Ehrenfest equations 
of motion for the electron density matrix and the mean value of the external spin (Landau equations): 
Due to the spin coupling they turn out to be nonlinear and, importantly, also coherences between 
electron states with different spin directions need to be considered. The electronic spin-polarized 
leads are implemented in form of a Lindblad-type dissipator in the infinite bias limit. We have 
solved this involved dynamical system numerically for various isotropic and anisotropic coupling 
schemes. For isotropic spin coupling and spin-polarized leads we study the effect of current-induced 
magnetization of the attached spin and compare this with a single quantum dot setup. We further 
demonstrate that an anisotropic coupling can lead to a rich variety of parametric oscillations in the 
average current reflecting the complicated interplay between the Larmor precession of the external 
spin and the dissipative coherent dynamics of the electron spin. 

PACS numbers: 73.23. Hk, 73.63.Kv, 75.76.+j, 85.75.-d 



I. INTRODUCTION 

Electronic quantum coherence in a few-orbital conductor 
[such as coupled semiconductor or molecular quantum 
dots (QDs)] which is subject to single-electron transfer—, 
can cause intriguing measurable effects such as transient 
current oscillations j 2 . complete current suppression^ and 
enhanced current fluctuations.— £ Such types of effects 
have been studied theoretically and experimentally in 
various QD geometries: in serially coupled dots^i 7 - in par- 
allel QDsj£ in triangular setups^ and in coupled orbitals 
in single molecules.— 

Of particular interest is the role of interactions of trans- 
ferred electrons with additional bosonic or fcrmionic 
degrees of freedom. In serial semiconductor QDs the 
coupling to a phononic bath can give rise to spon- 
taneous emission and absorption of phononS ) 11 ' 12 elec- 
tronic decoherence^^ or phonon replica in the transport 
characteristics^. Another type of coupling is provided 
by the hyperfine interaction with nuclear spins of the 
surrounding material^ - — This interaction enables spin- 
flip transfer between the electronic and the nuclear spin 
system, which resolves the current in the spin-blockade 
regime of serially coupled QDs. Moreover, in this regime 
it can cause coherent current oscillations with a period of 
the order of seconds^ which is believed due to dynamical 
polarization of the nuclear spins addressed theoretically 
in RefA 

In transport through molecular QDs 2 ^ the coupling to 
two types of degrees of freedom are of particular rel- 
evance: molecule vibrations 2 ^ and local magnetic mo- 
ments in single molecular magnetsi 22 ' 23 The latter es- 
tablishes the field of molecular spintronics 2 ^— and 
has attracted a large amount of theoretical work,— ~— 
which has focused on a single orbital as current-carrying 



state. For example, Bode et a/.— have derived Landau- 
Lifshitz- Gilbert- type equations of motion for the molec- 
ular spin based on a noncquilibrium Born-Oppcnhcimer- 
approximation, which assumes that the time scale of the 
dynamics of the local magnetic moment is much larger 
than the dwell time of the electrons. In Ref.— we have 
treated the dynamics of the average electron spin semi- 
classically in the infinite bias limit. 

In the present work, we will extend our method in Ref.— 
for the description of an electronic setup with two seri- 
ally coupled QDs, which interact with the same magnetic 
moment. We will show that in contrast to the single-QD 
case, the derived equations of motion for the electronic 
part will possess a much more involved structure due to 
the inclusion of coherences. The resulting coupling of the 
semiclassical dynamics of the attached magnetic moment 
with the spin dynamics of the non-equilibrium electrons 
gives rise to complex transient dynamics, such as chaos, 
and peculiar steady-state dependencies on the initial con- 
ditions. Furthermore, we find the phenomenon of para- 
metric resonance in the current oscillations. We compare 
the current-induced switching behavior of the attached 
magnetic moment in our coupled setup with the single- 
QD setup in Refs^i^. 

The remainder of the paper is organized as follows: 
Sec. HT1 contains the model with the Hamiltonian (|H AJ) . 
with a discussion of the level of description in Sec. Ill Bt 
with the semiclassical approximation and derivation of 
the Ehrenfest equations of motion pi Cp . and with the 
final equations of motion containing the transport dis- 
sipator in Sec. Ill Dl In Sec. IIIII the results will be dis- 
cussed: steady-state currents (|III A[) . isotropic coupling 
and current-induced switching of the large spin (|III Bl) . 
and anisotropic coupling (jlll Dj) . Finally, we will conclude 
in Sec. ITVl 
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II. MODEL 



Hamiltonian 



T, 



s R 



FIG. 1. (Color online) Setup of two QDs mutually tunnel- 
coupled with strength T c . The electron spins Si interact with 
a common large spin J via Xf . The QDs are attached to ferro- 
magnetic electronic reservoirs with spin-dependent tunneling 
rates Via-. 
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= 2i'* S 



(2) 



with 



the usual commutation relations 
= 5ij2 Sj and 5*, , Sj — ±<5y S • . 

In contrast to most of other related studies, here we 
also allow for anisotropic coupling between electronic and 
large spin: Xf ^ Af for a ^ f3. 

The lead electrons are assumed to be nonintcracting; the 
corresponding operator cj (q p<t ) creates (annihilates) 
an electron of momentum p and spin a in the Zth lead, 
while tip is the spin-independent coupling strength be- 
tween Ith lead and Ith dot. 



We consider a system of two serially coupled QDs with 
one orbital level each (see Fig. Q]). The system is subject 
to an external magnetic field B in the z— direction, which 
splits the QD spin levels. Moreover, the system is coupled 
to electronic leads, as well as to a large spin J with length 

J given by J |m, J) = J{J + 1) |m, J). The total system 
Hamiltonian reads as follows: 

H = H D QD + %3 + Pleads + H int (1) 

HDQD = £l die + 4j dja + B ^ S • 

i—L,R i7^j\c i 

Hj = B J , Hint — $i ^ 

a— x,y,z 
i 

Hlcads = 51 £ip<r cjpa CIP° + Yl + h ' C ' ' 

Here %dqd describes the double-QD system (DQD), 
with the tunnel coupling T c between the left and the 
right dot, and the energies of the orbital levels El,£r 7 
respectively. Note that our definition of B comprises the 
Bohr magneton and g factors. For the sake of simplicity 
we assume identical g factors for electron and large spin 
throughout the following. The Coulomb repulsion be- 
tween excess electrons within the QD system is assumed 
to be much larger than all other system and bath ener- 
gies to constrain their maximal number to one (i.e., the 
system is operated in the Coulomb blockade regime). 

The operators dl a (di a ) describe the cre- 
ation(annihilation) of an electron with spin a = "f, J, on 
the ith dot = L,R); the occupation operator is 

defined by h ia = d ia d i(J . The relationship between these 
electronic operators and the ath component of the spin 
operator Si is given by 



B. Level of description and time scale separation 

The microscopic dynamics of the system described by ([1]) 
is involved and contains a lot of information, which we 
are not interested in. Therefore we will choose a level 
of description, where only the dynamics of single-particle 
observables such as components of the average electron 
and large spin and the reduced density matrix of electrons 
dwelling in the dots will be considered. 
A rigorous derivation of the corresponding equations of 
motion in a controlled manner can be performed, e.g., by 
projective methods^ or by Kcldysh-Green functions^ 
The price one has to pay for projecting out some "irrele- 
vant" (microscopic) information is the occurence of non- 
Markovian terms and residual forces in the equations of 
motion. One typical way to retain Markovian (time local) 
dynamics may be the separation of time scales present in 
the considered problem. 

In this work we avoid this derivation path and gain the 
dynamics of the relevant observables by a simple semi- 
classical approximation instead (see Sec. Ill Cp . We ar- 
gue phenomenologically that the time scale of the large 
spin precession is much slower than the electron dwell 
time; i.e., electron spin fluctuations do not affect the large 
spin dynamics and vice versa. Those electron spin fluc- 
tuations arise from the tunnel coupling to the electron 
reservoirs, which we are going to treat in the standard 
Born-Markov approximation. 

In our Refj^ we have used a technique based on Laplace 
transform (see Appendix A of that reference) to combine 
the semiclassical Ehrenfest equations of motion for the 
spin observables with the reduced density matrix for the 
electrons. The same final equations of motions can be 
simply obtained by adding the terms based on the Lind- 
blad master equation for dot-lead coupling to the spin 
dynamics as we have proven for the single-QD setup in 
Refi2£. Here, we take advantage of this observation and 
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employ this simplified method to set up our final equa- 
tions of motion (see Sec. Ill Dj ). 



Semiclassical approximation and Ehrenfest 
equations of motion 



tion on the level of the interaction Hamiltonian: We sub- 



stitute Si =(S i )+5J and J = (J ) 
interaction Hamiltonian Ti- m t and obtain 



5 J into the 



(4) 



In order to investigate the dynamics of the coupled spin 
system we employ an equations of motion (EOM) tech- 
nique for the expectation values of the involved spin op- 
erators. In general, the EOM for the expectation value 
of an arbitrary operator O reads (h = 1) 



dt 



(O) 



0,H 



. .80. 



Due to the interaction between electrons and large spin 
in (fTJ) an infinite series of coupled EOMs of higher-order 
spin correlators will be obtained. In order to truncate 
this hierarchy we implement a semiclassical approxima- 



where the product of the spin fluctuators 5 S i 5 J is 
neglected, which is justified for J ^ 1 and Xf/B <C 1 
(Va, i). Thus, we essentially treat J as a classical object. 
Note, that since the Hamiltonian (fTJ) docs not contain 
any interaction of the large spin with the lead electrons 
or with an additional bath its length J will be conserved 

" ~2 



(3) on the microscopic level ( 



J ,H 



0). Nevertheless, 



later on we will introduce a way of damping in the large 
spin EOMs, which is not microscopically motivated. 
Using the semiclassical interaction Hamiltonian Q yields 
the following correlators for the electronic part of EOM 
(for the sake of clarity we omit the time dependencies) 



^ (4 4> = \ [ (<W - <W) { B + A * ( f )) + - <W) { B + % (J z )) + 2 - ] (4 *i°>) 



+ xT c [{dld jol )-{dld 3(T ,) 

+ \ [A? (f) i (*ot - 6„±) Af (?)] (4 4v) - \ [Af (f) i (5^i <W) A} (J*)] (4 4*') 



(5) 



where » = L(i?) for i = i?(L), ct =| (f) for <r =| (1). 
The EOM for the expectation values of the large spin's 
degree of freedom take the form of a Bloch equation com- 
pleted by the electronic back-action and the phenomeno- 
logical damping: 



d7< J » f > = 



Be z 



E 



$>(*) 



x (J) (t)—rj(J) (t), 
(6) 



with (S,) (t) ee J2 a Xf (S t ) (t)e a . 

So far we have derived the Ehrenfest EOM for the closed 
coupled spin system. In order to pursue a microscopic 
derivation of spin-dependent rates for the transport be- 
tween the electronic contacts and the system we could 
follow the method proposed in Refi2£. This technique 
requires the direct computation of all dot-lead corrcla- 



) and (Cipa>d jtT ), ro- 



tors appearing in Eq. ([5]): (d icr c 

spectively. However, the number of relevant dynamical 

variables in the DQD {d i<7 dj a ') (t) is far too large for 
the analytical derivation along the lines of Appendix A 
in Rcf>2£.. Instead we make use of the fact that for the 
single-QD the same equations of motion can be obtained 



alternatively by replacing the last term of the right-hand 
side of ([5]) with terms originating from the Lindblad-type 
master equation. Even though there is no proof of equiv- 
alence we believe that those aproaches lead to the same 
results in the case of DQD. The corresponding derivation 
will be provided in the following section. 



D. Transport master equation 



Here, we will combine the Ehrenfest EOM ([5]) with the 
quantum master equation of the DQD that is not cou- 
pled to the large spin. Specifically, we use the Lindblad- 
type master equation for the system density matrix p, 
which is derived by means of the standard Born-Markov- 
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approximation^ in the infinite bias limit: 



with p{t) ee ((p Q ) (t), P Ht),p±(t),Ut),Ut)) T and 



— p(t) = -i[H DQD ,p(t)} 



2 T La({dLad i L (7 J( t )} 



2 d La p (t) d La 



t b<t ({<4<t d-Ra- , } - 2 /5 (i) 



(7) 



;t 



Thc matrix elements of p are (d ia dj a ' ) . The tunnel rates 
are given by Ti a = 27r|i;| 2 pi a with the spin-dependent 
density of states pi a in the Zth lead, which is approxi- 
mated to be energy-independent. The asymmetry in the 
density of states will be parametrized by the degree of 
spin polarization (see also 4 ^) 



Pi 



_ Pit ~ Pij 
Pit + Pli 



(8) 



with pi 6 [—1,1]; pi = corresponds to a nonmagnetic 
lead and pi = ± 1 describes a spin up/down-polarized 
half- metallic ferromagnetic contact, respectively. The 
tunneling rates then read = |r;(l + pi) 1 Ti± = 

^r,(i- P ,), andr; = r> + iv 

Now we define the density matrix p in vector form with 
the following subvectors (omitting the time dependence) : 



p u = ({pI)Apr),{p1r)Aprl)) t 

= {(h La ) , (h R(J ) , (d La . d Rr7 ) , (d Ra d La ) 
L=((s1),01),0 + r),(Sr}) T 



(9) 



tu = ((d\t dm) . (d R i d L t) , (d m d Li ) , (d Li d Rt ) 
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where 044 is the (4,4)-zcro-matrix, O4 is the 4- 
dimensional zero-vector, V = (Tl<j, 0, 0, 0) T , 

r" ee (o,r fi(T ,o,o) T . 

Together with the modified Bloch equations © for the 
large spin (J) this forms an involved set of nonlinear dif- 
ferential equations, which describe the coupled dynamics 
of the non-equilibrium QD system and the large spin. 
Note that due to the internal tunnel coupling coherences 
between different spin directions in different QDs (£ e ,£u) 
also need to be considered in the EOM. The matrix on 
the right-hand side of (fT2|) provides the coupling struc- 
ture; the detailed definition of the 4x4 block matrices 
L' 7 , A, B, . . . F can be found in Appendix [3] In partic- 
ular, for vanishing coupling Xf = the p part decou- 
ples completely from the £-part and from the large spin 
(J). The resulting master equation for p describes the 
unidirectional transport of electrons through two spin 
channels ! 41 ' 42 Due to the strong Coulomb blockade the 
electron transfer in those channels is correlated and for 
finite coupling Xf they mutually couple due to spin-flip 
transfer between electron and large spin. 
Since we consider the QD system in the infinite bias limit 
the average time-dependent electron current through the 
DQD can be calculated by the product of occupations 
and tunnel rates of the right QD 



(I)(t) = T m (n R )(t)+T Ri {n R )(t). (13) 
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III. ANALYSIS AND NUMERICAL RESULTS 



A. Steady-state currents 



and the vacuum operator p = 1 — N with the total elec- 
tron number operator N = J2 la ni a ; its EOM is obtained 
by (J7J and reads 

^ (p Q ) (t) = -r L ( Po ) (t) + r«t (pi) (t) + r Ri (p R ) (t). 

(10) 

Eventually, the combination of ([5]), ©, and (|TU|) results 
in the EOM of the electronic part, which we provide in 
a compact vector form 

f t p(t) = M^{J x ),(J y ),{J z ))p(t) (11) 



Before we will start with the discussions of the dynam- 
ics of large and electronic spin, we consider the steady- 
state currents in the uncoupled and coupled case. With 
the interaction to the large spin turned off (i.e., A = 0) 
the steady-state current can be obtained analytically and 
reads 

(/) = 4r r 2 (r £i + r Lt )r JU iy 

e ' 4aT 2 + & i + 6 t ' 1 ' 

where a ee 2(T L1 -T R i + T Li T R ^) + r fii r flt , b a = 

^Lo^rs (r^CT 2 + 4e 2 ), and e = el — £r denotes the level 
detuning between the left and right dot levels in the sys- 
tem. With the definition of the tunneling rates and the 
assumption of symmetric tunnel couplings = T R = T 
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FIG. 2. (Color online) Current-induced magnetization of 
the large spin. (Top) Components of the electron spin {S L ) 
(blue, solid), (S R ) (blue, dashed) (S L ) (red, solid), (S V L ) (red, 
dashed), (§1) (black, solid), (S R ) (black, dashed). (Middle) 
Components of the large spin: (J ) (blue), (J v ) (red), (J ) 
(black). (Bottom) Electronic current (7) (t) vs. time. In the 
long-term limit one observes (J ) = (J y ) = 0, (J ) = J and 
the current vanishes. (Inset) Evolution of the large spin as 
trajectory in the Bloch sphere, where the red circle depicts 
the initial position. Note that the onset time depends on the 
initial orientation of the large spin, the more the large spin 
deviates from alignment with the magnetic field, the faster 
the switching process sets in. Parameters: pl = —pa = 1, 
B/F = 0.1, T c /F = 0.5, A/r = 1, e = 0, 7.//T = 0, large spin 
length J = 10. 



this can be rewritten as 



4T 2 (l-p|) 



(Jl = 

eT 4T C 2 [5 - p R {Ap L + p R )\ + 8e 2 (l - PlP r) + cT 2 ' 

(15) 



with c = (1 - 
tions p = pl 



PlPr)(1 — -P/f ) /2- For symmetric polariza- 
= Pr we further obtain 



eT 



AT?. 



20T 2 + 8e 2 + (l+p 2 )r 2 /2' 



(16) 



which becomes maximal for nonmagnetic leads (p = 0). 
We need to note that this current expression is only valid 
for incomplete polarization (|p| < 1) since in the polar- 
ized case the steady-state version of the master equation 
(|11[) is not well defined: In particular, the relevant part 
of the coefficient matrix (Liouvillian superoperator) M. 
(fT2")l reduces to the left upper NxN submatrix (detached 
by lines) with N = 9. For \p\ < 1 its rank is (N — 1) and a 
unique steady state follows. In contrast for |p| = 1 when 



one spin channel is completely switched off the rank re- 
duces and more than one steady state is obtained, which 
is clearly unphysical. The steady-state master equation is 
over-determined. In order to obtain the steady state for 
complete polarization we have to make use of the spin- 
polarized master equation considered in Rcfs. 5,42 where 
N = 5. The corresponding current reads 



eT 



4T 2 



12T 2 



r 2 ' 



(17) 



which significantly differs from (fTB]) and leads to a discon- 
tinuity of the current when \p\ approaches one. There- 
fore, the complete symmetric lead polarization provides 
a singular case in our description, which has to be either 
excluded or treated with caution (see below). 
Coupled spins - For arbitrary spin couplings we have not 
been able to derive an analytical expression for the cur- 
rent due to the high dimensionality of the EOM. How- 
ever, when the large spin components (J x ) and (J v ) ap- 
proach zero in the long-term limit, we again obtain the 
steady-state current (fT4|) . This either occurs when the 
large spin is damped (jj > 0) or for isotropic spin cou- 
pling (Af = A, Vi, a) as will be discussed in the next 
section. 

As in the uncoupled case complete symmetric contact 
polarization (|p| = 1) causes pecularities, such as a de- 
pendency of the stationary current on initial conditions 
or multistable current behavior (see Sec. IIII D]) . We note, 
that this is also caused by the above mentioned failure of 
the electronic master equation and not by the scmiclassi- 
cal approximation of the spin interaction. However, the 
spin-polarized master equation can not be used here. 



B. Isotropic coupling and current-induced 
magnetization of the large spin 

In this section we consider isotropic coupling Xf = A 
(Vi,a) between the electron spin and the large spin. We 
will discuss the dynamical behavior of both spins and the 
current for various lead polarizations pl and pr. 
Reverse lead polarizations. - For the sake of clarity we 
start our discussions with complete reverse polarization 
Pl = —pr = ±1. An electron transfer through the QD 
system only takes place when each spin-up (spin-down) 
electron entering the left QD is able to down-flip (up- flip) 
its spin state, respectively. According to the microscopic 

spin coupling Hint = ^^2j(^ §j J +5 Sj J + Sj J ) 
this electron spin-up (spin-down) flip is accompanied 
with the decrement (increment) of the large spin mag- 
netic number by one: m=Fl (J |m, J) = m |m, J) with 
|m| < J). Once the minimal (maximal) number =pj is 
reached, the large spin is completely aligned with the 
z— direction and the current vanishes. This process is in- 
dependent of the external magnetic field and can be con- 
sidered as current-induced switching^ of the attached 
local magnetic moment. This effect is also known in 
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(a) (b) 

FIG. 3. (Color online) Contour plots of the stationary currents vs. pl and pr. (a) Coupled QDs with parameters e — 0, 
T c /Y = 0.5, B/T = 0.1. (b) Single-QD with the same magnetic field B. The vertical arrows denote switching directions in 
areas divided by the straight dashed lines. For pl — Pr — 0, i.e., unpolarized contacts, and on the dashed straight lines no 
switching takes place. Note that for the DQD the transition between switching up and down depends on the product AJ: For 
AJ/r < 1/2 and XJ/Y > 2 the transition takes place along the line pr = pl; otherwise the transition line is located in the 
region bounded by pr ~ 2.6pi for AJ/r = 1 (dashed white line). Note that those transitions are not linear with respect to pl 
as shown, e.g., for XJ/Y = 1.25 (solid white line). 



magnetic layers described by a Landau-Lifshitz-Gilbcrt 
equation for the layer magnetization^ In this approach 
the current-induced magnetization has been interpreted 
as spin-transfer torqu o 43 ' 44 acting on the magnetic mo- 
ment. Unfortunately, for our DQD we are not able to 
identify analogous terms in Eq. [6] due to the complex 
coupling structure with the electronic part. 
Figure IIII Al shows the dynamics of the current-induced 
switching of (J) for our DQD, in the case of pl = —Pr = 
±1 and no damping of the large spin 7,7 = 0. In the tran- 
sient regime, damped coherent oscillations of the electron 
between the QDs with frequency 2T C (e = 0) are visible 
in the current and in the electron spin evolution. For 
the sake of clear demonstration we have chosen the ini- 
tial large spin to be reversely aligned to its final direc- 
tion. According to our above explanation the current 
then is peaked around the time when the derivative of 
(J ) (t) is maximal, in the isotropic case when ( J ) = 0, 
as the spin-spin interaction and thus the spin-transfer be- 
comes maximal there. In the long-term limit the current 
is given by (fT4|) . If one assumes a finite damping 7j > 
the z— component of the large spin will approach zero as 
well. 

The inset of Fig. IIII Al contains the large spin evolu- 
tion in the Bloch sphere - in addition to the large spin 
magnetization from ( J ) = — J to + J a precession with 
frequency B occurs. The time for the magnetization re- 
versal (switching time), however, does not depend on the 
external magnetic field. 



Since the coupling parameter A provides the number of 
spin-flips per unit time and T the number of electrons 
entering the left QD per unit time the switching time is 
determined by their ratio: It diverges when X/Y — > and 
saturates for X/Y ~ 1 as long as the tunnel coupling T c 
between the QDs is on the order of Y. This provides a 
lower bound for the switching time, whereas for T c /Y 3> 1 
and <C 1 it increases. 

Aside from the magnetization time scale, we observe that 
the onset of the switching process depends on the initial 
orientation of the large spin. In particular, if it is nearly 
aligned with magnetic field ((J } w J) the spin-flip rate 
becomes significantly diminished, since the correspond- 
ing coupling is proportional to A((J ) + (J V )) [sec (gj)]. 
This leads to a slowing down of the magnetization onset; 
in the limit (J } = J the large spin keeps his initial 
orientation. 

Arbitrary lead polarizations. - So far we have addressed 
complete reverse lead polarizations. However, switching 
also takes place for arbitrary pl and pn as shown in 
ITHAT a). The direction of the spin torque is schemati- 
cally depicted by vertical arrows depending on pi and 
Pr together with the steady-state current (QjJ as con- 
tour plot. For pi, < pr (pl > Pr) the large spin will 
switch to the spin-down (spin-up) direction, respectively. 
Note that this holds only for parameters XJ/Y < | or 
AJ/r > 2, respectively. If | < XJ/Y < 2 the transi- 
tion between up- and down-switching lies in the region 
bounded by pr — pl and pr ~ 2.6pL for XJ/Y = 1. In 
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FIG. 4. (Color online) Time-dependent currents in single- 
QD and DQD for different contact polarizations pl and pn 
during the process of magnetization reversal. The currents 
possess a maximum when ( J ) = 0. With decreasing \pl\ = 
\Pr\ the time for magnetization switching increases. In DQD 
the switching time is always larger than for the single-QD. 
Parameters: \/V = 1, T c /T = 0.5, B/T = 0.1 



this area the transition lines are not linear with respect 
to pl, as shown for XJ/T = 1.25 in Fig. IIII A7 a). 
The switching time is also affected by the choice of lead 
polarizations. Given that \pi\ — \pr\ — \p\, the switching 
time decreases with decreasing |p|, as depicted in FigJU 
Given a fixed left lead polarization \pl\ < 1, while pr is 
variable, the switching time has a lower bound for the 
minimum of —pl/pr- On the other hand the spin trans- 
fer is increased as the left lead polarization is increased. 
Thus, the switching time is minimized for complete re- 
verse lead polarization pl = —Pr = ±1. 
In the nonmagnetic case pl = Pr. = [origins in 
Figs. IIII Af a) and (b)] no switching takes place and the 
electron spin state ends up in a completely mixed state. 
When both QDs are presumed to be initially unoccu- 
pied the electronic spin state is completely mixed for the 
entire evolution. Consequently, the large spin evolves in- 
dependently of the electron spins. This enables an effec- 
tive description of the electron spin dynamics presented 
in Appendix |B] The situation, however, changes when 
the QDs are initially occupied. During the decay of the 
electronic state into a complete mixture the large spin 
evolution becomes affected and possesses transient oscil- 
lations before it runs into the back-action free precession 
for 7j = 0. The frequency of the transient oscillations de- 
pend on the parameter X/T while their duration is gov- 
erned by r. It follows that in the undamped case the 
large spin steady state depends on the initial QD occu- 
pation. Isotropic coupling with difference between left and 
right site - For a coupling scheme where the components 
of one electronic spin are isotropic, i.e., Xf = A; (Va) 
but the coupling differs between the two electronic sites, 
i.e, Xl 7^ Xr, we observe in principle the same current- 
induced magnetization behavior as before. However, the 
difference A A = Xl — Xr has an impact on the specific 
behavior, for instance, on the switching times. The more 
the couplings differ, i.e., for increasing |AA|, the longer 
the switching takes. Given that the coupling strengths 
are on the order of magnitude of the tunneling rates, the 
sign of AA is of subordinate importance. 



C. Comparison with single-QD 

In the following we will compare the phenomenon of 
current-induced switching in a singlc-QD and in the DQD 
in the same semiclassical description. The total Hamil- 
tonian for the single-QD (SQD) setup reads 

H = HSQD + Hi + Hint + Pleads (18) 

%sqd = X! £ dl d a + B S 

<r=t4 

Pleads = £ lP° 4pa CIP° + Yl U P £ U ^ + h - C - ' 

l,P,(7 lp <7 

where %j and Hint have been defined in (p}. Using the 
EOM technique introduced in Sec. Ill CI and the Lindblad 
master equation for the SQD setup 



at 



p(t) = -i[HsQn,p(t)] 

-\H Tl ° (K ^ >Ht)}-i£p (t) d a ) 

a 

+ r Ra ({(£ du , p (t)} - 2 d a p (t) dl 



(19) 

we obtain the following EOM for the QD occupations, 
the electron spin, and the large spin (r = Tl = Tr, 
omitting the time dependence): 



A (n°) = (a- (f) (s v ) - xy (f) (S X )) (S^ <W) 

-(TR a +T L a){n a )+T La 



dt 

d 
~t 



(S) = [(J) + Be z ] x (S)-T(S) 

-\r{p L +p R ){N) e z + ^Tp L e z 
[J) = [ (S) + Be z ] x (J) 



(20) 



with (S) = j: a X a (S )e a , (J) = J2 a X " (J ) e « an d 
(N) = (rvf) + (ni). The anisotropic version of these equa- 
tions (X x = X z = X and X y = 0) has been already studied 
in Refi2£, but with a different focus. For isotropic cou- 
pling Fig. IIII Af b) presents the steady-state current and 
the final direction of the large spin depicted as vertical 
arrows in dependence on the lead polarizations pl and 
Pr. The transition between the two directions occurs for 
Pl = Pr-, which has been also obtained in Ref. — . There 
the magnetization change was discussed in terms of the 
sign of the spin-transfer torque, which is determined by 
sgn[rL|F/?-|- — r£-fTflJ. In contrast to the switching in 
DQD [see Fig. IIII Af a)] the transition does not depend 
on A or J. We further note that since the electrons are 
assumed to be nonintcracting in the singlc-QD the cur- 
rent is symmetric with respect to an exchange of pl and 
Pr. This does not hold for the DQD. 
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The comparison of the current evolutions for both se- 
tups (as shown in Fig. [3J reveals that the SQD switching 
occurs always faster for the same set of parameters. In 
other words the coherent electron transfer between the 
QDs slows down the switching. 




80 160 240 320 

Ft 



FIG. 5. (Color online) Time-dependent currents in the case 
of undamped large spin for different initial system state: (a) 
empty DQD , (b) n Lt (0) = 1 , (c) n Li (0) = 1, (d) n m (0) = 1 
and (e) n^(0) = 1. The dominating oscillation frequencies 
are multiples of 2B for all currents (parametric resonance; see 
text), but we observe phase shifts and different magnitudes. 
Parameters: p L = p R = -0.9, B/F = 0.1, T c /F = 0.4, A£ = 
= = F, X R = 0. 



D. Anisotropic coupling 

Large spin switching. - Given the investigations of the 
spin magnetization reversal in the isotropic case we probe 
for anisotropic coupling schemes that yield spin switch- 
ing as well. In principle we observe that the switching 
process takes place as long as the a;-and y-components of 
the spins are isotropically coupled (Af w \\, Vz) and the 
leads polarized suitably. Taking into account the micro- 
scopic spin interaction [see Sec. IIII Bj this becomes clear 
since the A 2 terms are not directly contributing in the 
spin-transfer process. However, the z— couplings affect 
the speed of magnetization reversal: The magnetization 
process, e.g., can be slowed down by increasing the z— 
couplings on both sites: \ Z L — Afj = A 2 > A. Contrari- 
wise, decreasing A 2 steps up the switching. Even if the 
z-couplings are different switching occurs. The corre- 
sponding currents display a peak at the time when the 
derivative of (J ) (<) is maximal. 

Non-magnetic leads pl = PR = 0. - Basically the same 
as for the isotropic coupling holds. In addition we find 
that the transient phase and, thus, the polar angle for the 
free precession depend crucially on the initial electronic 
states. 

Undamped large spin 7,7 = 0. - We obtain various 
regimes with respect to lead polarization and spin cou- 
pling that exhibit limit cycles in the current, similar to 
Ref^. The Fourier spectra of those current time evolu- 
tions exhibit peaks at well defined frequencies, confirming 
the periodicity of the system's evolution. We note that 
in contrast to the single-QD^i the initial system occupa- 
tion plays an important role in the evolution of both the 



DQD and the large spin, which is illustrated in Fig. [5] 
A similar observation has been made in RefJ^. For a 
more detailed analysis the fix points of the EOM have 
to be determined. Unfortunately, due to the size and 
complexity of the EOM, we did not manage to find them 
analytically. Even their numerical computation provides 
a complicated task. 

Parametric resonance. - Studying the rich dynamics 
provided by the anisotropy, we can observe parametric 
resonance in the current oscillations, particularly. The 
evolutions exhibit current oscillations at multiples of the 
doubled Larmor frequency of the large spin, i.e., 2B, as 
can be seen in Figs. EJa)-[6jc). Although the periodic 
oscillations are in general very complicated, due to the 
nonlinearity in the equations, the dominating frequency 
of the oscillations of the x and y components of the elec- 
tron spins, respectively, is B. The z component and con- 
sequently the current, on the other hand, oscillates with 
multiples of 2B, which we reveal in the Fourier spectrum 
[Fig. Eld)]. The same phenomenon has been observed in 
the single-QD, so that we can employ its model to pro- 
vide an explanation. We assume anisotropic spin-spin 
coupling, i.e., X x = A 2 = A and A y = and magnetic 
leads (pl,Pr 0)- Given that A J <C 1 we can make 
use of the back-action free EOM for the large spin (|B.2I) . 
It follows from Eqs. ([20) that (tf) and (S V ) oscillate 
with frequency B. We further recognize that the time 
derivative of {n 7 ) couples to the product (S V ) (J }. Inte- 
grating this product of sinusoidal oscillations both with 
frequency B leads us to the current and (S ) oscillating 
with frequency 2B^. 

Given that the large spin is damped by the rate 77 the 
evolutions for different initial DQD occupations match in 
the long-time limit. Consequently, the periodic or nonpe- 
riodic oscillations are only present during the transients 
(as shown in the panels of Fig. [5]). We observe different 
transient scenarios ranging from quasiperiodic [(a)-(d)] 
to chaotic motion [(e)-(h)]. In all cases the steady-state 
current is given by (|14|) . 

IV. CONCLUSIONS 

In this work we have studied electron transport through 
a DQD setup coupled to electronic contacts . The spin 
of the excess electrons in the DQD interacts with a large 
spin and an external magnetic field is applied. We use 
semiclassical Ehrcnfcst EOM together with a quantum 
master equation technique. This method works well if 
one assumes that the dwell time of the electrons is much 
smaller than the average large spin precession period. 
We have found that the coupled dynamics of the large 
spin and the electron spins as well as the current through 
the DQD structure strongly depend on the polarization 
of the electronic leads and on the isotropy of the spin- 
spin coupling. Particularly, if the coupling between elec- 
tron spin and large spin is isotropic with respect to x 
and y while the electronic leads arc reversely polarized, a 
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setup with strongly polarized contacts (pl = Pr = —0.9 and 
3 ) and different couplings between spins: (Left) \1 = Ar = \ Z R — Y, A^, = and 

Ayr 
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0.7. (a) and (e) Components of the electron 



FIG. 6. (Color online) Different transient currents for 
T c /r = 0.4, B/Y = 0.1, e = 0, jj/Y = 10 
(Right) AJ/T = 0.2, A*/r = 4.2, A£/r = 1.3, AS/T = 2.9 
spin (5*i) (blue, solid), (S R ) (blue, dashed) (S L ) (red, solid), (S L ) (red, dashed), (S^) (black, solid), {S R ) (black, dashed), 
(b) and (f) Components of the large spin: (J ) (blue), (j v ) (red), (J ) (black), (c) and (g) Electronic current (7) (i) vs. time. 
The steady-state current is given by Eq. (|14[l. (d) and (h) In the current spectrum of (c) the dominating frequencies are even 
multiples (2,4,6. . . ) of the large spin Larmor frequency B, which refers to the phenomenon of parametric resonance (see text). 
The current spectrum of (g) is broad and noisy, which we attribute to chaotic behavior (see also Ref.— ). For both cases holds: 
As the large spin is damped out, the current oscillations vanish, since no more spin-flips occur. 



current-induced magnetization process of the large spin 
is obtained. The speed of this switching is always lower 
than for a single-QD setup and depends, most impor- 
tantly, on the lead polarization. 

On the other hand, a complete anisotropic coupling leads 
to a rich variety of different dynamical scenarios, which 
for an undamped large spin depend on the initial state. In 
particular, we have found that the electronic system may 
show self-sustained oscillations, either quasi-periodic or 
even chaotic. Introducing a large spin damping renders 
these dynamics to be transient. A remarkable feature 
of the undamped or damped dynamics is the occurence 
of parametric resonance, observable as frequency dou- 
bling in the current. Furthermore, we analytically have 
derived the stationary currents for the spin-dependent 
unidirectional single-electron transport, when no interac- 
tions with the external spin are present. It turns out that 
this current also applies for the case of current-polarized 
or damped large spin. 

Experimental realizations for the large spin could be ei- 
ther the spin of magnetic impurities in semiconductor 
QDs or a net magnetic moment in single molecules (single 
molecular magnets) . Whether our mean- field approach is 
applicable for these examples, where the spin typically is 
not very large, needs further investigations. However, one 
possible realization for the large spin, that might justify 
the semi-classical treatment is the hyperfine interaction 



with an ensemble of nuclear spins, where the number of 
spins is reasonably high4^ 
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Appendix A: Block-matrices and parameters 

Here, we provide the coupling matrices used in the EOM 
for the electronic part (fTT]) 



{J ) 



/ o 





-iT c 


iT c 
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iT c 


-iT c 
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iT c 


-Zf- + ie a 







V iT c 







Fro- 
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ie a j 



{(f) , (f)) 



(-K A£ 
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-A* A+ 
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B ( (J ) , (J ) ) = 



((f) , (f) 



E 



f-iT c 
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D y(J ) 

f((j 2 ) 
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v° 


-A£ 


A fi 













-iT c 
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- A R 

° y 



-iT e 



-iT c 








= Diag^ie L , -i£ L ,i£R 
= Diag I ie' + - 



2 ' 



2 



, — lei ^ i ie ^— i 

2 + 2 2 



At 



2 ; 



with the parameters: 



1 



e^e + -{S^-5 ai ){\l-y k ){J ) 
e' ± ^±e + B + l -{\l-Y R )(f) , 



(A.l) 



it = B + Xf (J ) ,(i = L, R) 



A 



± _ 



Af(J )±iA?(/> 



(A?)* = -A? 



Appendix B: Large-spin dynamics with vanishing 
electronic back-action 



Without the electronic spin the EOM of the large spin 
© or in flU) read 



— (J) = B x (J) - 7J (J) 



(B.l) 



which are readily solved by 



'cos(-Bi) -sin(St) 0^ 
(J)(t) =e~ ljt I sin(Bt) cos(Bt) I (J) Q . (B.2) 







1, 



This provides the damped Larmor precession of the large 
spin around the external magnetic held axis with fre- 
quency B. Inserting Eqs. (|B.2|) into the electronic EOMs 
leads to a set of linear first order differential equations 
with time-periodic coefficients. In two dimensions it cor- 
responds to the Mathieu-Hill equation, which describes 
the phenomenon of parametric resonance. 
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